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ABSTRACT 

Matter in the Universe is arranged in a cosmic web, with a filament of matter typically connecting each neighbouring galaxy 
pair, separated by tens of millions of light-years. A quadrupolar pattern of the spin field around filaments is known to influence 
the spins of galaxies and haloes near them, but it remains unknown whether filaments themselves spin. Here, we measure 
dark-matter velocities around filaments in cosmological simulations, finding that matter generally rotates around them, much 
faster than around a randomly located axis. It also exhibits some coherence along the filament. The net rotational component is 
comparable to, and often dominant over, the known quadrupolar flow. The evidence of net rotations revises previous emphasis 
on a quadrupolar spin field around filaments. The full picture of rotation in the cosmic web is more complicated and multiscale 
than a network of spinning filamentary rods, but we argue that filament rotation is substantial enough to be an essential part 
of the picture. It is likely that the longest coherently rotating objects in the Universe are filaments. Also, we speculate that this 


rotation could provide a mechanism to generate or amplify intergalactic magnetic fields in filaments. 
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1 INTRODUCTION 


The cosmic web of bubble-like voids, separated by walls, filaments 
and haloes, describes the matter distribution of our Universe on scales 
much larger than galaxies. It is predicted by the standard model 
of cosmology (Klypin & Shandarin 1983; Bond et al. 1996) and 
observed in galaxy surveys (de Lapparent et al. 1986; Colless et al. 
2001; Zehavi et al. 2011). Intergalactic filaments are the skeletons of 
the cosmic web, generally connecting pairs of neighbouring galaxies. 
They feed matter into nodes of the cosmic web, and halo and galaxy 
spin is known to correlate with orientation with respect to filaments; 
small haloes' spin vectors tend to align with the axis of filaments 
they inhabit, while large haloes' vectors tend to be perpendicular to 
the axis (Aragón-Calvo et al. 2007b; Hahn et al. 2007; Paz et al. 
2008; Tempel & Libeskind 2013; Dubois et al. 2014; Codis et al. 
2015; Wang & Kang 2017; Ganeshaiah Veena et al. 2018; Wang 
et al. 2018; Wang & Kang 2018; Kraljic et al. 2020). This finding 
for small haloes is one motivation for our study: small haloes deeply 
embedded in a spinning filament might tend to spin along with it. 
The current understanding of spin in the cosmos begins with the 
tidal torque theory (TTT, Hoyle 1949; Peebles 1969), a formalism 
that has seen rather continuous conceptual study (e.g. White 1984; 
Porciani et al. 2002a,b; Schafer 2009; López et al. 2019; Motloch 
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et al. 2021; Neyrinck et al. 2020; López et al. 2021). In the TTT 
approximation, a collapsing object is approximated by an aspherical 
blob of matter with a moment of inertia that generally does not align 
with its tidal tensor. This misalignment causes the blob to torque 
up as it collapses, an idea that applies equally well in 2D (i.e. to a 
filament cross-section) as 3D (i.e. to a halo, the usual context for 
these ideas). 


The same tidal field that is involved in spin generation can also be 
used to understand the formation of both filaments, and the haloes 
that typically form at their endpoints. At these haloes (large density 
peaks), the tidal field generally prescribes inward flow in all direc- 
tions. Along density ridges that join these nearby peaks, the tidal field 
generally prescribes filament formation, i.e. expansion along the axis 
between the peaks, but (nearly cylindrical) collapse along the two di- 
rections perpendicular to it (van de Weygaert & Bertschinger 1996; 
Bond et al. 1996; van de Weygaert & Bond 2008). 


At the same time as a filament collapses, the tidal field typically 
produces a quadrupolar pattern in the vorticity and ‘spin field’ around 
it (Pichon & Bernardeau 1999; Codis et al. 2012, 2015). The spin 
field (r x v, with r being the vector between the filament axis and 
a particle, and v its velocity) around an ideal filament points along 
the filament axis; it alternates sign across four quadrants, separated 
by the primary and secondary axes of collapse. Importantly, the 
total angular momentum (which we also call ‘spin’), summing over 
the four quadrants, is usually ignored. The spin is precisely zero if 
integrating over a cylinder in the irrotational initial conditions, but 
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is generally nonzero at later times. The question we ask is whether 
filaments spin substantially in the late-time Universe. 

Filamentary rotation has been speculated about before. There has 
been work on filaments providing a 'swirling, rotating environment 
that ends up spinning up nearby haloes (Codis et al. 2012; Laigle 
et al. 2015). In an art-inspired toy ‘origami approximation,’ haloes 
spin if and only if filaments attached to them do (Neyrinck 2016b); 
this suggests that since galaxies generally spin, so do filaments. Also, 
haloes spin in 2D cosmological simulations; extruding them to 3D 
suggests that filaments spin substantially (Neyrinck et al. 2020). But 
a convincing measurement of how coherently realistic 3D filaments 
spin is still missing. 

This paper is organized as follows. In $2, we describe the Mil- 
lennium Simulation (MS) dataset and methods for our fiducial 3D 
measurement, using an observationally motivated filament definition. 
We use the MS because it is a standard dataset in the community, 
with a good balance between statistics and mass resolution. In $3, we 
give these measurements. In $4, as a guide to understanding these 
3D results, we show the dynamics in 2D, of *halo' (roughly, filament 
cross-section) spin in a 2D simulation. This 2D evolution can be 
tracked much more easily and visually than in 3D, and readers may 
even wish to read this section before looking in depth at $82 and 3. 

In $5, welook in further depth, and put our results in the context of 
a physically motivated filament definition involving (nearly) cylindri- 
cal collapse between haloes. For this analysis, we use two further 3D 
simulations, both with slightly attenuated small-scale modes, which 
clarifies the dynamics. One of them uses the initial conditions of the 
Illustris simulation, which will be helpful in our future studies look- 
ing at properties of gas and galaxies in filaments. The wide range of 
simulations used in our study all have similar results, indicating their 
robustness. 

In $6 we give some thoughts on a somewhat subjective issue, of 
why we consider the ‘longest coherently rotating objects in the Uni- 
verse' to be filaments. Finally, in $7, we conclude with a discussion 
of observational prospects and implications of our results. 


2 METHODS 


For our fiducial measurement, we use a simple filament definition, de- 
signed to make straightforward contact with observations: the straight 
line connecting pairs of dark-matter haloes above a mass threshold. 


2.1 The Millennium simulation 


The Millennium simulation (MS; Springel et al. 2005) is a dark- 
matter-only cosmological N-body simulation, using 2160? dark mat- 
ter particles to represent the matter distribution of a model Uni- 
verse in a cube of 500 Mpc/h on a side. Each dark-matter parti- 
cle has mass 8.6 x 10847! Mo. The simulation was run in a fidu- 
cial ACDM model with the following cosmological parameters: 
Qm = 0.25, Qp = 0.045, h = 0.73, Q4 = 0.75, n; = 1 and og = 0.9. 
For fast computation, we randomly downsampled to one-tenth the 
original particle density. 

Dark matter haloes were detected in the simulation with a Friends- 
of-Friends (FOF) (Davis et al. 1985) method, resulting in the same 
halo sample as in the Millennium database (Lemson & Virgo Con- 
sortium 2006). We define a halo's position as the location of its 
most-bound particle, involving a Spherical Overdensity approach. 
See Cai et al. (2017) for further details. We use haloes with mass 
greater than 10? 57! Mo to mimic a sample of Luminous Red Galax- 
ies (LRGs) (Parejko et al. 2013; Zhu et al. 2014), assuming that each 
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Figure 1. Illustration of coordinates for the halo-filament system. The line 
connecting the two haloes at ry and r» is defined as the filament axis. The 
filament region is defined to be 1 Mpc/h away from each halo, to isolate the 
filament signal from possible contamination from flows within each endpoint 
halo. The average angular momentum Jayg (see Eq. 1) for an individual 
filament defines the direction of the spin. 


halo hosts one galaxy. We choose this galaxy type because it is 
well-characterized observationally, the target of many large-scale- 
structure surveys (e.g. Dawson et al. 2013). We end up with 35,300 
haloes in the box. Our fiducial filament catalogue consists of the 
33,951 halo pairs separated by 6-10 Mpc/h. This distance range is 
chosen so that the halo centres are substantially farther apart than 
their own ~ 1 Mpc// radii, but are near enough to be typically con- 
nected by a coherent dark-matter filament (Colberg et al. 2005). For 
longer separations, there are many more pairs (3,732,470 pairs for the 
longest separation we considered). To keep the noise level similar for 
all separations, we analyse 30,000 randomly chosen pairs in each of 
the longer cases; adding more filaments changed results negligibly. 

What does a ‘coherent dark-matter filament’ mean? There is not a 
universally agreed-upon, absolute definition of this (Libeskind et al. 
2018), but our conception, shared by many, is a density ridge (e.g. 
Sousbie 2011) between two haloes. Density ridges are thought to 
exist even without haloes at both endpoints, but these would be 
difficult to detect observationally for modest-size filaments, so we 
only consider filaments with two halo endpoints. Generally, within 
high-enough density ridges, dark-matter streams have crossed, i.e. 
initially distant dark matter particles have come together in the same 
place, moving at different velocities, possible for such a collisionless 
fluid (e.g. Falck et al. 2012). Also, that gas has shocked (the baryonic 
version of stream crossing). We expect that a typical filament in our 
fiducial catalogue (defined as straight-line segments between LRGs 
in this length range) corresponds approximately to some density- 
ridge filament. This sort of definition has been used extensively in 
observations to detect observational signals related to filaments, such 
as from lensing, and the thermal Sunyaev-Zeldovich effect (Clampitt 
et al. 2016; Epps & Hudson 2017; de Graaff et al. 2019; Xia et al. 
2020; Tanimura et al. 2019). This definition is particularly useful 
when there may be no other reliable tracers of the structure between 
haloes available. For comparison, we also test longer filaments in the 
following ranges: 10-20, 20-40, and 40-60 Mpc/h. 


2.2 Filament definition and stacking 


As shown in Fig. 1, for each pair of haloes at positions ry and r2, we 
define the filament axis as the line connecting them, with the direction 
n= aa For every particle with position vector x, L = (x r1): f 
is then the projected distance from r4. 

We define the filament region to be the dark matter between the two 
haloes, and at least 1 Mpc// away from each halo along the filament 
axis to avoid confusion from velocites within each halo. The average 
virial radius r200 of our halo sample is 0.5 Mpc/h, ~ 1 Mpc/h in 
the 99th percentile. Halo regions are the matter within 1 Mpc/h from 
their centres along the filament axis. For each particle in the filament 
region, we calculate the projected angular momentum Jproj = (r x 


v) - fi. We divide particles into 15 equally spaced cylindrical shells 
around the filament axis, ranging from 0 to 6 Mpc/h. Here, we choose 
the spin direction A = +1 for a filament so that the average angular 
momentum over the cylindrical bins is positive. That is, that 


-2 y-1 
A D shells i (Zoiie k dogs) Ci N; 
Javg = —— n > 0 (1) 
=) > 
2 shells i 0j 


where N; is the number of particles in the ith cylindrical shell, c2 
is the variance of J over particles in the shell, and the subscript k 
labels each particle within the shell. To stabilise Jayg statistically, we 
inverse-variance-weight the sum over shells. 

Having defined coordinates, we measure filaments’ density and 
angular-momentum density fields in rectangular boxes along the axis, 
and then project these boxes either along or perpendicular to the axis. 
To cleanly combine results from different filaments, for each filament 
we rescale the distance along the filament r|. so that the two haloes 
are centred at +1, but keep the filament/halo region boundary fixed 
at 1 Mpc/h (Fig. 1). 


2.3 Measuring rotational velocity and angular momentum 
density 


We measure the rotational velocity of dark matter in cylindrical shells 
for the filament and halo regions. The average rotational velocity in 
the ith shell is defined as 
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where N;. y is the number of particles within the ith shell in the fth 
filament, and r, ; is the perpendicular distance to the ith shell. For 
the average angular velocity in haloes, we repeat above calculation 
but for particles that are within 1 Mpc// to the haloes. 

We use the term angular momentum density for a quantity con- 
ceptually defined as J = (1 + 6)vrotr1, i.e., a quantity that volume- 
integrates up to the total angular momentum of the filament (divided 
by the mean density). Note that this is not the same as the curl of the 
momentum-density field (1 + 6)v, which is more directly related to 
the vorticity. Explicitly, we measure J in the ith shell as 

2 filaments f 2 particles k J proj, ksi: f “A 
One =< ET TU = (3) 
2 filaments f Rsep;f Aip 
where Rsep;f is the length of the filament region, A; is the area of 
ith annulus, and p is the mean density in the simulation. 


3 RESULTS 


Fig. 2 shows the projected matter and momentum density p = (1 + 
ó)v, averaged over 33,951 filaments 6-10 Mpc/h long. (1+ ô = p/p 
is the matter density p relative to the mean matter density p, and v is 
the velocity.) In addition to the radial infall due to gravity, there is a 
clear rotational component around the filament axis. These two sum 
up to the spiral pattern shown in Panel A. In the next two lengthwise 
panels, we rescale the coordinate along each filament axis to put all 
endpoints at +1. (B) shows that, on average, the pattern is coherent 
across the filament. In (C), for each halo pair with 4 = —1, we flipped 
the sign of velocities instead of interchanging positions of the pair, 
before stacking them. This procedure aligns the sense of rotation but 
flips any infall velocity near the two haloes. When averaged over a 
large sample, this has effect of nearly nulling flows along the filament. 
The rotational velocities close to the two haloes are small; it seems 
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that for this sample of rather large haloes, rotations of their host dark 
matter haloes are not particularly aligned on average with that of 
the filaments. In the filament region (between the haloes, but with 
distance along the filament >1 Mpc/h from both — see Fig. 1), the 
average rotational velocity peaks at around 80 km/s near the filament 
axis. In (D), the rotation pattern has a peak at 1Mpc/h, and thus 
is neither like a solid body (with velocity increasing with distance), 
nor like a fluid vortex (with decreasing velocity). Also of note is 
that, in the halo region (outside the filament region), the velocity 
it is much smaller. Thus the endpoint haloes only slightly corotate 
with the filament. We find that for the innermost bin, the period of 
rotation is of order the age of the Universe, befitting the Universe's 
longest spinning objects. This Hubble timescale also accords with the 
magnitude of vorticity generated in caustics (Pichon & Bernardeau 
1999). The rotation frequency has an exponential distribution, going 
up to ten times the mean in our sample. Still, this timescale is long 
compared to the crossing time of a typical filament; it seems that just 
like in haloes, radial motions dominate rotational motions for the dark 
matter, i.e. dark-matter filaments are not typically rotation-supported. 

We repeated the measurement for filaments with the least and most 
massive 30% of halo pairs (judging by the average mass of the pair). 
Curiously, within r} ~ 1 Mpc/h, the rotation is highest for low-mass 
pairs, but this relationship flips outside this radius. It is tempting 
to think of this inner regime as *within the filament, but typical 
filaments could be curved on scales up to ~ 1 Mpc/h, so we caution 
against this conclusion. 

To test if the filament-rotation signal might arise from random 
velocity flows, we repeated the measurement around pairs of random 
positions with the same separations as the haloes. The average angular 
momentum density profile J (r , ) in the random filaments is positive, 
to be expected since it is on average aligned with Jays in the same way 
as in the real sample, but in the randoms it is much smaller. Panel E 
of Fig. 2 shows that J (r , ) from filaments is significantly higher than 
Zero, after subtracting the signal from random filaments. For a plot 
showing both reals and randoms without the subtraction, see Fig. 11, 
using a warm dark matter simulation, with similar results as in the 
MS. In summary, angular momenta of real filaments are much larger 
than those of randoms. 

We also investigated how the signal varies with filament length 
(Panel E). We find that the rotation decreases for longer filaments, 
becoming very small at 40-60 Mpc/h. This is expected, since pairs of 
haloes at this large separation rarely have nearly straight-line physical 
density ridges connecting them. 


3.1 Relation to quadrupolar flows 


This coherent filament rotation has not before been measured, to our 
knowledge. How is it related to the quadrupolar spin field mentioned 
above (Codis et al. 2015)? This finding is also related to the quadrupo- 
lar vorticity field expected around filaments at late times (Pichon & 
Bernardeau 1999; Laigle et al. 2015). (Recall, however, that we are 
showing the angular momentum density, rather than vorticity here.) 
It arises from the curl-free, gravitationally-sourced velocities of an 
elliptical density peak ($3.1); also see Fig. 7. The orientation of each 
filament cross-section is random in Fig. 2, averaging out any such 
pattern, so in Fig. 3 we orient the stack along the major axis (the 
largest eigenvector) of the projected matter distribution around each 
filament. The density concentration along the x-axis indicates a wall 
on average, within which the filament is embedded (left panel in 
Fig. 3). The overall spin is positive everywhere, with its amplitude 
varying with orientations (middle panel). 

Subtracting the isotropic signal (measured from the left panel of 
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Figure 2. Top panels: Average density and momentum density fields around filaments. Colors show the average density around filaments connecting pairs 
of haloes separated by 6-10 Mpc/h found in the MS. Fig. 1 is a sketch of the coordinate system. (A) Arrows show the momentum density p in the filament 
region, projected along the filament axis, excluding the halo regions. The longest vector corresponds to 1380 km/s. (B) Lengthwise view on one side of the 
filament axis, stacked over filaments (rescaling coordinates so each halo lies at +1), and projected. Here, arrows show v instead of p, for clarity at large 1+ ô. 
The longest vector corresponds to 264 km/s. (C) Same as (B), but with infall velocities around two haloes nulled (See text). Bottom panels: (D) Rotational 
velocity profile, and mass dependence. Average azimuthal velocity vrot as a function of distance to the filament axis r1, in the filament region and the halo 
regions as shown in panel B. Matter in the filament region (blue) exhibits stronger rotation around the filament axis than in halo regions (black). Dotted and 
dashdotted blue curves show the profiles for low- and high-mass subsamples, defined to contain filaments in the lowest and highest 30% in the mean mass of 
the two haloes. (E): Filament angular momentum density profile, J (r1 ). The upper panel compares J (r1) for different filament lengths in units of Mpc/h. 
Signals from random patches of the same length in the simulation have been subtracted for each case. The lower panel shows their corresponding error bars, 
which are the standard errors of the mean velocity of a filament in the filament stack. (F) Rotation frequency. The mean frequency at each distance, excluding 
any filaments that counter-rotate (compared to Jzyg) in that distance bin. (G): Rotation frequency distribution over filaments, in their innermost bins (blue). 
Excluding the 23% of filaments that counter-rotate (compared to Javg) in this bin, the mean frequency is 1.05, i.e. a rotation period of about the age of the 
Universe. An exponential distribution (orange) with this mean fits the distribution well. 


Fig 2) reveals a quadrupole pattern (right panel in Fig 3). Importantly, 
the average amplitude of the monopole, which was not considered in 
the past, exceeds or is comparable to the average quadrupole. This 
suggests that it should not be ignored when interpreting the align- 
ment of galaxy spins with their filaments. The white zero regions 
separating the quadrupole regions grow vertical and horizontal at 
large radius, as expected (Codis et al. 2015). But closer to the fila- 
ment, the quadrupole is ~ 45° rotated, perhaps from the matter flow 
overshooting the filament with some impact parameter. 

To illustrate a simple situation in which a quadrupolar pattern of 
angular momentum arises, we show a 2D elliptical density peak in 
Fig. 4. Following Eq. (12) in Mandelbaum et al. (2006): 


2A a|, B(a- 2)? 
(a — 2)2 ala — 4) 
with parameters œ = 0.8, ej, = 0.3, A = 1000 and B = ey,a/2 for 


$(r,0) = cos(20) | , (4) 
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illustration. The convergence x is then computed via 


I 2 
= -V’o 
K=5 (5) 


using Fourier transforms. 

We use the Zel’dovich approximation (Zeldovich 1970) to estimate 
the gravitationally-sourced velocity from this, taking the gradient of 
the inverse Laplacian of the 2D density field, i.e., v = V®. Note 
that these velocities are exact in the initial conditions, and should 
be good approximations on large scales at late times. Finally, we 
compute the angular momentum density using the modelled density 
and velocities. The result is shown in the right panel of Fig. 4; a pure 
quadrupole with no monopole. A typical quadrupole arises simply 
from the potential flow that an elliptical density peak produces. This 
has no monopole in the initial conditions, and the monopole that we 
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Figure 3. Panel A: Average density (colors) and momentum density fields (vectors) around MS (panels A-C) filaments, as in Fig. 2A, but stacked so that the wall 
in which a filament might be embedded will appear horizontally. The maximum length of vectors corresponds to 1760 (km/s). Panel B: Magnitude of angular 
velocity with the positive direction defined as into the paper. Panel C: Angular velocity with the isotropic component subtracted. The isotropic, monopole 
component dominates the quadrupole. Panel D: Monopole (Jo) and quadrupole (J2) amplitudes as a function of distance, averaged over 32 filaments of length 
6-10 Mpc/h from a rerun of the Illustris simulation. These followed the collapsed ridge of the filament, rather than the straight line between endpoint haloes. 
For the orange ‘filament-aligned’ curve, Jz was measured for each filament; for the green 'segment-aligned' curve, Jz was measured separately in each of 6 
segments along each filament, resulting in increased Jz amplitude in the inner bins. See $5 and Fig. 14 for further discussion. 
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Figure 4. An idealized elliptical density peak (Left). The gradients of its 
gravitational potential give rise to a quadrupolar pattern of angular momentum 
density (Right). 


see at later times comes from fluctuations away from this in halo 
shape, in the density field or in its Lagrangian boundary. 


3.2 Rotation signal from random samples 


Our measure of angular momentum is by construction positive, so 
its positive average is meaningless. To attach physical meaning, it is 
essential to know the amplitudes of the signal arising from a random 
scenario. For this, we define pairs of random points chosen to have the 
same lengths as the real filament samples. We repeat measurements 
as shown in the previous sub-section for the angular momentum 
density profiles for these random samples. For most of the cases we 
have investigated, the random signal is found to be sub-dominant over 
the signal from real filaments. The profiles shown in Fig. 2E have 
had the random signal subtracted off. For the longest filaments (40- 
60 Mpc/h), the residual J(r  ) is very small, indicating that the signal 
is very close to the random case. See Fig. 11 for a comparison of 
typical real and random signals, without subtracting off. In estimating 
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Figure 5. The cumulative distribution functions (CDF’s) of the average an- 
gular momentum Javg for random line segments and filaments. The green 
dashed line indicates the maximum difference between the CDF's, the D 
statistic used in a two-sample Kolmogorov-Smirnoff test. The black dotted 
lines indicate the thresholds corresponding to 68.3% (1-0) and 95.5% (2-07) 
probabilities that Jayg is not drawn from the Jayg distribution of random line 
segments. In the random sample, the fractions exceeding these thresholds are 
1 — 0.683 = 31.7% and 1 — 0.955 = 4.5%, to be compared respectively with 
79% and 26% of real filaments. 


the error bars shown, we neglected the small error on J(r , ) measured 
from the randoms. 

A standard, single summary statistic of the difference between 
such distributions uses the two-sample Kolmogorov-Smirnov test. 
We find a maximum vertical separation between the CDFs of D = 
0.47. For this sample of tens of thousands, the null hypothesis that 
the two samples are drawn from the same distribution is rejected 
with essentially certain confidence; the p-value characterizing the 
probability that the two distributions are the same is ~ 10-120, 

More tangibly, how might we answer the question of whether an 
individual filament is substantially rotating or not? We can decide by 
reference to the distribution of Jayg for the random sample, shown in 
Fig. 5. If we define a filament to be ‘rotating’ if its Jayg is over a 1-0 
cut (where 31.7% of the randoms have higher Jayg), 79% of the real 
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Figure 6. Filament-to-filament correlation matrix of the angular momentum 
density J (r1) in different radial bins. Each sample is a single filament with 
length 6-10 Mpc/h. The positive off-diagonal elements indicate that matter 
in different radial bins indeed generally co-rotate. This correlation decreases 
with increasing bin-separation, as expected. 


filaments are rotating. If we are more strict, and define ‘rotating’ with 
a 2-0 cut (where 4.55% of the randoms have higher Jaye), 26% of 
the real filaments are rotating. This 2-0, 26% fraction is our fiducial, 
simple response to the question ‘how often do filaments substantially 
rotate?' but it uses an admittedly arbitrary threshold. 


3.3 Rotation coherence with radius 


The direction of angular momentum for each filament is defined 
with an inverse-variance-weighting scheme as described above. The 
question remains how coherent the rotation is with radius. To check 
this, we compute the correlation matrix for the rotational velocities 
r X v at different radii. An example for the 33,951 filaments with 6- 
10 Mpc/h separation is shown in Fig. 6. We find that the correlation 
coefficients between adjacent radial bins are close to unity, indicating 
that the matter in adjacent cylindrical shells strongly tends to co- 
rotate. The coefficient decays with bin separation, as expected. 


4 2D SIMULATION 


In this section, we examine halo collapse and rotation in 2D, as a 
conceptual guide to filament collapse and rotation in 3D. This corre- 
spondence is obviously only approximate. Dynamics in 2D and 3D 
can qualitatively differ in many physical systems. But the relatively 
simple dynamics of the cosmic web on megaparsec and larger scales 
encourages us to think that looking at a simplified 2D setting here can 
help, as long as we keep in mind the caveats that we are neglecting 
variation along the filament, and motions out of the plane (such as 
filament bending, or rotation perpendicular to the axis). 

The monopole we find is not present in the primordial, irrotational 
velocity field of the simulation. Integrating a zero-curl field around 
a circle, the angular momentum vanishes. Where, then, does the 
rotation at a later time come from? It comes from the non-circular 
outskirts of a filament cross-section in the initial conditions, generally 
carrying angular momentum (Doroshkevich 1970). To understand the 
collapse of matter onto a cross-section of a filament, we additionally 
investigate a two-dimensional N-body simulation. 

This simulation was previously used by Neyrinck et al. (2020), with 
some analysis and animations (see https://youtu.be/7KjesL.. 
hP7c) relevant to the present work. The simulation has 1024? parti- 
cles, and a box size of 32 Mpc/h. The initial conditions were a 2D 
slice of particles from a 1024? simulation with a BBKS (Bardeen 
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et al. 1986) initial power spectrum, with modes of wavelength 
< 1 Mpc/h suppressed in amplitude, to simplify the physical problem 
and clarify structures like filaments. The (x-y) plane of particles was 
replicated along the z-axis to match the number of slices in the orig- 
inal particle grid. It was run using a version of GApGET2 (Springel 
2005), modified to calculate forces and update positions only in the 
x and y directions. Because of this replication along the z-axis, each 
particle represents a cylinder, interacting with other cylinders effec- 
tively with a 2D version of gravity. So, each of its haloes is literally 
an infinite filament. For the initial power spectrum and expansion 
history, we used a generic ACDM (Qy = 0.3, Qa = 0.7, og = 0.8, 
h — 0.7) set of cosmological parameters. 

We detect haloes in the 2D simulation using the oRIGAMI (Falck 
et al. 2012) algorithm. A particle is classified as a halo particle if, 
going from the initial to final conditions, it has ‘folded’ (crossed some 
other particle) along two (in 2D) initial orthogonal axes. We then join 
together groups of particles adjacent on the initial Lagrangian square 
grid to form haloes.! 

Initially, the angular momentum has a clear quadrupolar pattern in 
the centre (Panel A of Fig. 7), summing to zero around circles. At the 
outskirts, though, the boundary of the protofilament is not circular, 
giving a net angular momentum (Panel B) (Doroshkevich 1970). As 
this matter falls onto the filament, pulled by gravity, it carries this 
angular momentum to the filament axis. Collapsed particles appear 
randomly arranged, losing their initially obvious quadrupole (Panel 
C). Here, ‘collapsed’ particles may not be virialized; they simply 
experience multistreaming along two orthogonal axes. In Panel D, 
we quantify the evolution of the monopole for a stacked sample of 
183 such haloes in the simulation. As time passes, the averaged an- 
gular momentum is transported from the outskirts toward the centre. 
Collapsed particles contribute most to the angular momentum near 
the centre of the filament, whereas non-collapsed particles dominate 
the angular momentum at the outskirts. Near the filament centre, the 
ratio between the monopole and quadrupole is zero initially, but then 
increases with time as the monopole increases and the quadrupole 
diminishes. In summary, angular momentum around filaments orig- 
inates according to the tidal-torque theory (Peebles 1969; White 
1984): an asymmetric matter distribution is torqued up and is later 
transported gravitationally to the centre. As with dark-matter haloes 
(White 1984; Motloch et al. 2021), the angular momentum of a fil- 
ament cross-section seems to grow as the scale factor a?/2 until 
collapse, after which angular momentum is conserved in physical 
coordinates, generally retaining its direction (Neyrinck et al. 2020). 


5 ANALYSIS IN TERMS OF PHYSICAL FILAMENT 
PROPERTIES 


Here, we investigate filament spin in greater depth, particularly im- 
portant because of our empirical, rather than physical, filament defi- 


! Several haloes returned in this process were actually groups of haloes 
apparently distinct by eye, joined by small bridges of halo particles. The 
spurious fragmentation of filaments into tiny haloes in simulations with trun- 
cated initial power (Wang & White 2007) likely contributed to this. Applying 
a mathematical morphology erosion operator cut many of these bridges, re- 
turning haloes that generally corresponded to visual expectation. Erosion, 
used e.g. by Platen et al. (2007), shaves off cells within a specified distance 
(here, one pixel) of a boundary from all contiguous blobs. We computed the 
erosion by smoothing the Lagrangian ‘halo’ (1) and *not-halo' (0) field with 
a circular top-hat filter of radius one pixel; after smoothing, we classified as 
*halo' particles all corresponding pixels with a value 0.99. 
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Figure 7. Evolution of halo angular momentum in a 2D simulation. If we neglect variation along a 3D filament, this halo from a 2D simulation corresponds 
to a filament cross-section. Colorbars in Panels A-C apply to the angular momentum J. In Panels A & C, particles in the simulation that collapse to this halo by 
z = 0 are shown with bold colours; other particles are shown with muted colours. (A) At z = 10, the quadrupole pattern is clear both for the protohalo region 
and the surroundings. (B) The cumulative net angular momentum of particles, summed from the centre, shows that J originates at the outskirts of the region 
that later collapses. (C) The quadrupole randomizes during evolution from z = 10 to the current epoch (for clarity, J for the non-collapsed particles has been 
transformed to sinh"! (J/10)/2.) (D) The average angular momentum in radial bins from 183 such 2D haloes at three different epochs shown in the legend. 
We see the angular momentum being transported to the centre. Here rii; ax is the largest comoving distance from the area centroid to a particle in the initial 
conditions, and r is in comoving coordinates. Solid curves include all particles, while dashed curves include only particles that are collapsed at z = 0. The 
shaded regions represent the error on the mean. The amplitude is lower than in Fig. 2E; note though the visible noise from the smaller sample here. 


nition. First, in $5.1, we look at how these filaments defined simply 
as regions between rather large haloes correspond to what we expect 
visually. $5.2, we look at how the spin correlates with the density 
and arrangement of matter in cylindrical filters such as used in our 
filament definition. 


In the rest ofthe section, we use alternative simulations, practically 
because they are accessible for analysis with the oRIGAMI web-type 
classifier that we use, but also this broadens the range of cosmologi- 
cal settings of our results. In $5.3, we look at the degree of physical 
collapse in empirical filaments. $5.4 examines the amount and co- 
herence of spin along collapsed density ridges between haloes, rather 
than along straight lines. For this analysis, we use both the oRIGAMI 
and MMF?2 filament-finders to find these ‘collapsed density ridges.’ 


Examining rotating filaments in terms of the oRIGAMI and MMF2 
definitions, in addition to our fiducial definition, should sufficiently 
elucidate the spin's physical nature and definition-dependence for 
the purposes and scope of this paper. But there are several other 
well-motivated filament finders in common use that are worth men- 
tioning. A non-exhaustive discussion of them follows; for a more 
comprehensive comparison, see Libeskind et al. (2018). 


Several filament-finders work with the density field or the point 
process itself of galaxies or haloes; this is the most straightforward 
class of filament-finders to apply to observed galaxies. DisPERSE 
(Sousbie 2011) has already been used to study the vorticity and spin 
field around filaments (e.g. Pichon et al. 2011; Codis et al. 2015). 
It finds density ridges, conceptually-similarly to MMF2 (see also 
the original MMF, Aragón-Calvo et al. 2007a) and Nexus (Cautun 
et al. 2013). But DisPERSE works on the scale of individual particles, 
rather than on the field as smoothed on multiple scales, which can 
cause results to differ somewhat. The Bisous filament-finder (Tempel 
& Libeskind 2013) also works with galaxies as points, finding strings 
of points commonly appearing together in cylindrical filters; Bisous 
filaments are broadly similar to density ridges. 


Another, dynamical class of algorithms brings in further infor- 
mation beyond the density. oRIGAMI is this category, using ordering 


differences between particles’ Lagrangian (initial) coordinates and 
Eulerian (late-time) coordinates. Others use the tidal force field con- 
figuration (Hahn et al. 2007; Forero-Romero et al. 2009), or the 
velocity shear field (Hoffman et al. 2012). See Ganeshaiah Veena 
et al. (2018) for some discussion of how filaments defined by the 
density and velocity shear field relate to each other. 


5.1 Visual filamentarity between Millennium halo pairs 


Although the density stack in Fig. 2 of 33,951 empirical filaments 
(defined as regions between two large haloes) looks quite like an 
idealized filament, we wondered how they individually correspond to 
visual expectation. Colberg et al. (2005) found by visually inspecting 
halo pairs in this length range that they are joined by a recognizable 
“straight or warped’ filament 86% of the time; however, they used 
haloes with mass threshold 10!4Mo /h; with our smaller 10M o/A 
cut, we expect the fraction of recognizable filaments in our sample 
to be somewhat smaller. 

Here we show flyarounds of the density fields around four example 
halo pairs from our sample, spanning a range of filamentarity: 1; 2; 
3; 4. Unsurprisingly, looking at several examples, we saw many cases 
where the pair of haloes are connected by a matter bridge close to a 
straight line, but also many cases otherwise. Nonetheless, it is clear 
that on average, there is a matter bridge in our sample (Fig. 2), and 
they rotate substantially. 


5.2 Dependence of spin on matter arrangement in a cylinder 


As one way to elucidate the physical nature of the filaments in our 
fiducial MS sample, we look at a filament's rotation speed depends 
on its average density. We do so by splitting our filament sample 
according to the mass density averaged within the cylindrical filament 
region, ô. As shown in Fig. 8, the rotational velocity vrot increases 
with ó. This is expected from the gravitational origin of the rotational 
velocity, since filaments with more mass have brought matter from 
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Figure 8. The dependence of the angular momentum per particle on the 
density within cylindrical ‘filament region’ filters, for six radial bins (indicated 
by titles above each plot). The density-rotation relation for random cylinders 
(one point per cylinder) is shown in blue, and for our sample of halo-bracketed 
filaments in orange. Dashed lines show best-fitting linear models. 
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Figure 9. The rotational velocity per particle as a function of distance to 
the filament axis. The blue line shows the result from pairs of haloes that are 
separated by 6-10 Mpc/h. The orange line shows the result from random pairs 
of points with this range of separation, but with a halo of at least 10? Mo /h in 
the middle.The pink line shows the expected value based on the mass-rotation 
relations found in Fig. 8. 


larger distances, with a longer lever arm in their angular momentum 
sum. In Fig. 9, we show that by using the best-fitting relation between 
Vrot Versus ô, the rotational velocity is accurately recovered. 
Curiously, though, Fig. 8 shows a density dependence even for 
randomly selected line segments. Is this all we are seeing in our 
results, that our filaments are just particularly dense regions using 
our cylindrical filter? We investigate the effect of the shape of the mass 
arrangement in Fig. 9. We select random ‘filaments’ between pairs of 
points centred on random haloes. Typically, this matter distribution 
would be far from filamentary. In the figure we find that the rotational 
velocity of the non-filamentary matter distribution is much weaker 
than the case of filaments. This suggests that it is important to have 
an actual filamentary matter distribution to give substantial rotation. 
In another test, we look at a projected cross-section of random 
cylinders in the top 10% of the density range. In Fig. 10, both the 
selected filaments and random line segments are on average cylinder- 
like structures, as in Fig. 2A. Whether we select cylinders that lie 
between a pair of haloes, or simply mass within that cylindrical filter, 
we tend to pick up filamentary regions. In this sense, our definition of 
filaments is rather general, and we expect the results to hold for a wide 
range of filament definitions that detect elongated mass distributions. 
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Figure 10. The average matter-density cross section of the top 10% densest 
random line segments, as measured in the inner 1 Mpc/h. 


5.3 Relation to a collapsed-filament definition 


Here we study an aspect of our filament definition, exploring the con- 
tribution of the spin from high-density regions that have undergone 
cylindrical collapse. In an idealized picture, the spin would come 
from such physically collapsed filaments, but with our empirical def- 
inition, we expect that instead, some halo pairs do not have collapsed 
density ridges between them. 

There are various possible degrees of ‘collapse.’ An idealized 
filament would have experienced cylindrical collapse, along two di- 
rections perpendicular to its axis (its particles classified as ‘filament’ 
by or1GAMI, which we use for this test). Realistically, most would 
contain haloes that have collapsed along all three axes as well, classi- 
fied as ‘halo’ by oRiGAMI. We count these as part of the filament. Such 
fragmentation would be entirely expected for even an idealized fila- 
ment in the MS, since small-scale fluctuations typically drive small 
haloes to form in such a filament. We also wish to test to what degree 
a filament sample like ours contains wall or even void regions. 

We could not apply this analysis to the MS because ORIGAMI can- 
not analyse a simulation with glass initial conditions (like the MS). 
Instead, we use a warm-dark-matter (WDM) dark-matter simula- 
tion from Yang et al. (2015), which assumed ACDM cosmological 
parameters slightly different from the MS. Smoothing the initial con- 
ditions produces a more visually prominent cosmic web (see Fig. 3 
of Neyrinck 2012). The 100 Mpc/h, 5123-particle simulation we use 
had initial conditions smoothed at a = 0.1 Mpc/h (but with a gradual 
kernel, extending a factor of ~ 6 larger in scale). 

We identified the haloes used to define filaments in the 3D WDM 
simulation in a new way, similar to that used for the 2D simulation: 
we first tagged particles with oRIGAMI (Falck et al. 2012). We joined 
together ‘halo’ particles that were adjacent on the Lagrangian grid. 
We first applied an erosion operator with a radius of one grid spacing, 
to cut some visually spurious links, and to remove tiny spurious 
haloes within filaments (Wang & White 2007). We used this method 
because the outer (‘splashback’) caustics of haloes in WDM are 
more visually evident than in CDM, and the oniGAMI method worked 
well to identify these halo boundaries (Neyrinck 2012). This fully 
Lagrangian method may be promising more broadly as a halo-finder, 
but is especially useful for a WDM simulation with a clarified cosmic 
web. The haloes we used had at least 1000 particles, resulting in a 
sample of 1825 haloes. Because of the substantial differences in the 
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Figure 11. Angular momentum density profile J (r, ) contributed from par- 
ticles with different origami tags in the WDM simulation, for filaments 6- 
10Mpc/h long (left) and for a random sample with the same length (right). 
The tags indicate the number of axes of collapse; i.e., O (void), 1 (wall), 2 (fil- 
ament — according to ORIGAMI, not our straight-line definition), and 3 (halo). 
J (r1) increases as the degree of collapse increases. Error bars indicate the 
error on the mean. 


halo finder, mass estimate, and WDM cosmology, we do not treat 
this halo sample as equivalent to our MS sample, but these results 
are still relevant. 

In Fig. 11, we see that for filaments in the WDM simulation (left), 
the signal is dominated by halo particles, with decreasing contribu- 
tions from filament and wall, and totally uncollapsed void particles. 
Note that the ‘halo’ particles are not in the haloes at the end of the fil- 
ament, but are in (likely smaller) haloes between them. Also, because 
of spurious fragmentation in WDM filaments (Wang & White 2007), 
some filament particles may have been classified as halo particles. 
So, we consider both ‘halo’ and ‘filament’ particles to be filament 
matter that has collapsed. Another reason to consider ORIGAMI ‘halo’ 
particles to be collapsed along the filament is for near-invariance to 
a small-scale power cut-off. If this were a CDM rather than WDM 
simulation (e.g. in the MS), the filament would have fragmented into 
haloes. At right, we do the same measurement, except using random 
filament endpoints instead of haloes. In that case, the contribution 
from each type is much smaller. 


5.4 Sensitivity to filament definition, and lengthwise rotation 
coherence 


Here, we take this analysis further, tracking filament spin closely 
along collapsed density ridges between haloes, almost never straight. 

Pichon et al. (2011) made a related measurement, in a paper fo- 
cusing on how filaments contribute to halo and galaxy angular mo- 
mentum. In their Appendix A, they show a correlation function of 
angular momentum vectors along DisPerSE filaments, finding that 
angular momentum along filaments is highly correlated on scales up 
to 15 Mpc/h. 


5.4.1 Method 


For this measurement, we analysed filaments identified in a dark- 
matter-only re-run of the Illustris (Nelson et al. 2015) simulation, 
assuming cosmological parameters from the Nine-Year Wilkinson 
Microwave Anisotropy Probe analysis (Hinshaw et al. 2013). We 
did not expect qualitative differences in our results from this slight 
change in cosmological parameters, but it is good to check. In our 
re-run of Illustris, we smoothed its initial conditions with a 2 Mpc/h 
sharp-k filter. We could have used this Illustris-rerun above in $5.3, 
but we used both, to explore results in a variety of simulations and 
halo-finding methods (there, orRIGAM1, applied directly to a WDM 
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simulation, and here, using large haloes identified with FOF in the 
unsmoothed Illustris simulation). 

To use ORIGAMI in this context, it was necessary to infer a set of 
grid (instead of glass) initial conditions from Illustris-3-Dark. We 
formed a regular 5123 grid by interpolating the velocity fields from 
the first publicly available snapshot (z = 46.77), using a Delaunay- 
based interpolation scheme (van de Weygaert & Schaap 2009). We 
then recovered grid initial conditions using the inverse Zel'dovich 
approximation to recover the Illustris initial conditions. We checked 
that by eye, the large-scale cosmic-web structure in our remapped 
runs and the original Illustris-3-Dark correspond. 

We generated an MIP ensemble (Aragon-Calvo 2016) of 32 re- 
alizations from the smoothed initial conditions, allowing negligible 
discreteness noise. The initial conditions of the 32 realizations share 
the same Fourier modes above a 2 Mpc/h cutoff scale, but are pop- 
ulated with different modes on smaller scales. Using particles from 
all realizations of this correlated ensemble allowed us to compute 
velocity fields with negligible discreteness noise for each filament, 
even with a naive cloud-in-cell velocity estimation method. 

We applied two different web classifiers to identify the collapsed 
density ridge: MMF2 (Multiscale Morphology Filter 2, based on the 
density Hessian, applied at multiple scales) and orIGAMI (a dynam- 
ical definition using directions of particle crossings). There were 86 
pairs of haloes of mass at least 10'3Mo/h that satisfied our separa- 
tion cut of 6-10 Mpc/h. We found that 32 of these were separated 
by clear bridges of MMF2 filament-classified voxels; all of these had 
clear bridges of orIGAMI particles. There were also typically small 
haloes (and particles oRIGAMI-classified as ‘halo’) between the end- 
points; we consider any haloes between the endpoints to be part of 
the filament, and include their particles in the analysis. 

We grouped particles into 7 quantiles along the ‘filament region’ 
(Fig. 1) of each filament, measuring the spin within each segment. 
We positioned the quantiles to give (almost) the same number of 
particles in each segment, rather than an equal-length division. 

We traced the ridge by linearly interpolating between nodes, one 
node per quantile. We positioned each node in x and y (if the z-axis 
is along the straight line separating endpoints) by finding the median 
x and y particle coordinates in that quantile (neglecting particles 
that have distance larger than 3 times the standard deviation in that 
bin). A ‘segment’ occupies the space between adjacent nodes; there 
are 6 segments between the 7 nodes of a filament. In each segment, 
we measured the angular momentum around an axis connecting the 
nodes that bracket it. An example of the segmentation is shown in 
Fig. 12. 

We then measured 2D rotational velocity profiles for each filament 
segment, setting the positive direction to align with the filament’s 
rotation within the inner 1 Mpc/h, averaged among all segments. To 
achieve negligible particle discreteness noise, we stacked particles 
over 32 MIP realizations. For each of the 2D rotational velocity 
maps, vg(r, 0), we measure the amplitude of the quadrupolar pattern 
at different phases ¢ at distance bin r; via 


2n 
O(ri, 9) = f T vo(r, 8) cos(20 + $)rdrd0, (6) 
0 reri 


and we define the phase of the quadrupolar pattern to be the value $1 
of 9 that maximizes the quadrupole ampliude Q, measured from pix- 
els within 1 Mpc/h of the filament axis. At each radius, we set J2 (ri) 
at radius bin r; to be the amplitude of the quadrupole with phase 9 
measured in the rotational velocity field vg(r;, 0), including pixels 
within the radius bin r;. Note that even though Jo(r < 1Mpc/h), 
averaged over the inner 1Mpc/h, is positive, J2(r) can go negative 
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Figure 12. Segmentation of a filament in our re-run of the Illustris simulation. 
Particles tagged by orIGAMI as filament and halo particles are coloured blue. 
The thick black curve shows the 6 segments of this filament. A fly-around 
animation of this example is available at https: //qx211.github.io/img/ 
Filament 65.mp4. 


at larger radius, if the phase of the quadrupole at that radius differs 
sufficiently. 


5.4.2 Results 


An example filament is shown in Fig. 13; a 2D rotational velocity 
map is shown for six segments. The quadrupolar (and more general 
anisotropic) pattern is indeed visually evident in this case. From the 
clock plot, we can also see the rotation of the quadrupole among 
the six segments, i.e. the direction of the quadrupole rotates from 
segment to segment along the filament. In Fig. 14, we can see that 
the amplitudes ofthe monopole and quadrupole in each segment vary, 
i.e. Jo can be larger or smaller than J2. But even though visually, what 
stands out is the anisotropic pattern, the monopole is dominant at low 
radius, in most segments. Once averaged over the entire filament, 
Jo is generally dominant over J2, which we show in Fig. 3D, for 
easy comparison with the Millennium results. This monopole can be 
compared to the MS filament velocity profile in Fig. 2D, which is 
larger; this is counter-intuitive if we expect the collapsed ridge (which 
we trace in Illustris) to have a greater speed. The small sample size 
(only 32 filaments) makes this a weak statement, however. 

In Fig. 15, we show average flows onto filament collapsed density 
ridges (at center) over this 32-filament sample, rotating the signal 
in each of the 6 segments to align as in Fig. 3. The results are 
qualitatively the same as in that figure, but the horizontal ‘wall’ 
in the density field is hardly evident, because of the much smaller 
sample. The velocity field displays similarity to Fig. 3C, though. 

We show the velocity field of a ~ 4Mpc/h-long filament from the 
Illustris simulation in Fig. 16. An animation of advected particles in 
another view is linked in that caption. What is only a hint of a helical 
loop in this frame is a much clearer helical flow in the animation. 
The velocity field in this animation is an average over DTFE (van 
de Weygaert & Schaap 2009) velocity fields from 32 Illustris MIP 
realizations, corresponding to an initial power truncation at 2 Mpc/h 
(explained above). The velocity field was further high-pass-filtered 
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with a 0.5 Mpc/h Gaussian, to remove large-scale flows. We used 
the Unity engine for the advection and rendering. 


5.4.3 Summary 


In this subsection, we used another standard reference simulation (Il- 
lustris, with many different realizations of small-scale modes), and 
two different filament definitions, going beyond our fiducial straight- 
line definition. By studying individual filaments in depth, we saw that 
the amplitude of the monopole and quadrupole along the filament 
are generally comparable, with the phase of the quadrupole shift- 
ing. When averaging over the entire filament, because of this phase 
rotation, the averaged quadrupole is typically suppressed compared 
to its value in single segments. So a straight-line filament definition 
may overestimate the dominance of the monopole over quadrupole, 
but qualitatively, the results are the same: the spin-field monopole is 
generally of comparable or higher magnitude than the quadrupole. 


6 CANDIDATES FOR THE LONGEST ROTATING 
OBJECTS IN THE UNIVERSE 


Here we clarify and elaborate on our statement that the longest sub- 
stantially rotating objects in the Universe are likely filaments. We 
word this carefully; we do not mean that all filaments are substan- 
tially rotating, and are longer than all other substantially rotating 
objects in the Universe. Instead, we mean that the few longest rotat- 
ing objects are likely to be filaments. 

An ‘object’ here is a cosmic-web component, a coherent virialized 
or collapsed object (halo, filament, or wall), or a region without 
collapse (void). For definiteness, we assume a dynamical dark-matter- 
sheet folding definition of these objects, as in ORIGAMI (Falck et al. 
2012), but some other definitions broadly agree, e.g. the MMF2 
(Aragon-Calvo & Yang 2014), and Nexus (Cautun et al. 2013); see 
Libeskind et al. (2018). In the oRIGAMI stream-crossing definition, a 
supercluster, an overdensity in the density field smoothed on at least 
tens of Mpc/h, would likely contain many virialized haloes, voids, 
collapsed filaments, and collapsed walls. 

A filament typically originates from an oblate, disk-like object in 
the initial conditions (Zeldovich 1970; Shandarin & Zeldovich 1989; 
Falck et al. 2012; Lovell et al. 2014; Hidding et al. 2016; Feldbrugge 
et al. 2018). The oblate object turns prolate, stretching along the axis 
that becomes the filament axis, and compressing toward it. As with 
a halo, any small angular momentum in the collapsing protofilament 
leads to faster rotation, like a figure skater pulling in their arms. 

But there is no equivalent process in the formation of a void or 
wall that is likely to speed up rotational motions. Walls do collapse 
along one axis, and can contain rotating patches (haloes and filaments 
within them). However, they are unlikely to rotate coherently (e.g. as 
a disk), since they generally expand faster than the cosmic mean in 
the two directions perpendicular to the collapse axis (Aragon-Calvo 
et al. 2011), evolving similarly to voids in a 2D Universe. This serves 
to suppress rather than enhance rotational motions. 

Voids are even less likely to rotate coherently. They are “cosmic 
magnifiers’ (Aragon-Calvo & Szalay 2013), their interiors expanding 
faster than the cosmic mean, resulting in vanishing vorticity (e.g. 
Hahn et al. 2015), and negligible large-scale rotation expected. 

As for haloes, they are well-known to spin, but are of maximum 
radius (e.g. as measured from the Millennium haloes) ~ 2 Mpc/h, 
and so do not compete in size with filaments. 

Another, mere suggestion for filaments as the longest objects is 
that in the toy, origami approximation (Neyrinck 2016a,b), filaments, 
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Figure 13. Left: Average rotational velocity field, in km/s for each of the six segments of a filament. Right: The average rotational velocity profile from combining 
particles in the six segments together. The clock plot in the lower-right panel shows how the quadruople's phase changes across segments. 12 o'clock is defined 
as the phase estimated from the entire filament. 
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Figure 14. Monopole Jo(r) and quadrupole Jz (r) from the velocity (in km/s) maps shown in Fig. 13. Blue line thicknesses and opacities increase with segment 
number. The phase giving the maximum amplitude for J is set in the innermost 1 Mpc/h, and then is used to measure J) at all radii. It goes negative where the 
quadrupole’s phase is sufficiently rotated from the inner region. In the middle panel, this phase is set by the vg averaged over all segments. In the right panel, 
the phase is determined segment-wise, i.e. to maximize Jz in the inner 1 Mpc/h of each segment. This is why the right panel generally has higher amplitude 
than the middle. 
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Figure 15. Similar to Fig. 3 A and C, but for a much smaller sample of 
32 Illustris filaments. Left: the stacked density, with arrows representing 
velocities. The maximum length of a vector corresponds to 1176 km/s. Vector 
lengths £ were transformed using £ = 50 sinh"! (v/50), with v in km/s; e.g. 
a vector representing 50 km/s is 4.4 times shorter than the maximum. Right: 
the azimuthal component of the angular momentum. In both stacks, each of 
the 6 segments along each filament has been rotated so that the PCA major 
axis of the density distribution aligns with the horizontal axis. Deep blue 
patches at right correspond to large clockwise velocity arrows at left. 


Figure 16. A view of particles (golden comets) advected along streamlines 
of the final snapshot of the velocity field in our smoothed re-run of the Illus- 
tris simulation (starting the animation after some advection has already hap- 
pened). We highly recommend watching our animation showing another view 
of the particles being advected, at https: //www.youtube.com/watch?v- 
hi-a-htHAxY 


which can be arbitrarily long, are the largest rotating components; 
walls and voids do not rotate. Voids do not rotate by construction in 
this model, but the non-rotation of walls is not trivial. 

The situation in the real Universe is more ambiguous, of course. 
While spin along filaments has substantial coherence, this coherence 
degrades with distance. So the longest filaments observed may have 
only slightly aligned spin at both endpoints, so the most extreme 
objects would likely be open to interpretation. 

Also, patches exist that straddle the clear categories of wall, fil- 
ament, and halo. Boundaries between objects (or objects and sub- 
objects) can be ambiguous, as well. While we have argued that walls 
and voids are unlikely to rotate substantially, it would be interesting 
and relevant to measure the degree of wall and void rotation in sim- 
ulations, although their complex geometries would make this a more 
subtle measurement than for filaments. 
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7 DISCUSSION 


We find that cosmic filaments generally carry some spin, which has 
some coherence along them. Within the filament, the dark matter has 
substantial random and radial motion, though; one should not think 
of filaments as nearly uniformly rotating cylinders. The velocity field 
in a filament consists of overlapping streams of dark matter going 
in different directions, with a net rotation after averaging over all 
streams. One might hesitate to call this rotation, but in fact this is 
the same situation as with haloes, which are said to rotate, despite 
even more severe multistreaming. We also should clarify that there 
is typically much substructure along each filament; no individual 
filament looks as idealized as the stack in Fig. 2. 


As with halo spin, we advocate thinking of filament spin as a 
fundamental property, even though that spin can be small compared 
to random motions, and a practical spin measurement for a filament 
(as with a halo) depends on details and definitions. But we must also 
acknowledge that even though it is generally a dynamically younger 
object, the velocity field of a filament can be more complex than of a 
halo. A halo is readily approximated by a single velocity and a spin, 
while in principle the spin along a filament can vary with position 
along it (typically mildly, as we have shown), and the filament’s shape 
could distort, something that we have not measured here. 

During the review process of this paper, we learned that another 
team (Wang et al. 2021) recently claimed detection of a sample of 
rotating filaments outlined by galaxies in the Sloan Digital Sky Sur- 
vey, prompted by our present work. Still, we give our own thoughts 
about detecting filament rotation in the observed Universe here. The 
most straightforward way to detect rotation in filaments may be to 
identify thick filaments of galaxies between clusters; one could look 
for rotation in the galaxies and gas between them such as using ra- 
dio surveys to locate turbulent gas in intra-cluster bridges (Brunetti 
& Vazza 2020). For a rotating filament nearly in the plane of the 
sky, on average, one side will move away from, and the other to- 
wards, the observer, as in Fig. 2B, giving a redshift and blueshift, 
respectively. That is, redshift-space distortions would tend to shear, 
or tilt, the filament along the axis, such that the side moving toward 
the observer is closer. This might be detected by stacking galaxy 
redshift-space distances along filament axes, and looking for dipolar 
asymmetries in each stack. Assessing this method’s feasibility may 
require a full mock observation, since it would sensitively depend on 
the galaxy sample, and the redshift-space dispersion within the fila- 
ment. It would be important to rule out any mechanisms that might 
mimic a dipolar asymmetry in galaxy positions, such as filaments 
tilted even in real space. This would be a substantial investigation, 
which we leave to future work. 


Going beyond galaxies, any ionized gas around the filament, likely 
co-rotating with the dark matter and galaxies, should scatter photons 
from the cosmic microwave background (CMB) in opposite direc- 
tions in the two sides, by the kinematic Sunyaev—Zeldovich (kSZ) 
effect. This would cause a dipole extending along the filament axes, 
imprinted on the CMB temperature map. It might be measured in 
a manner similar to how the kSZ effect has been detected around 
galaxies and galaxy clusters (Planck Collaboration et al. 2016; Boc- 
quet et al. 2019). Again, assessing this detection’s feasibility would 
be a substantial study, depending sensitively on both the amount of 
ionized gas rotation, and observational details, so we leave this to 
future work. 

Even purely theoretically, it will be interesting to study how gas 
rotates within filaments. As in haloes, the collisionless dark matter 
likely has substantial radial velocities in a filament; a typical position 
near the axis would have many streams going both in and out. Gas, on 


the other hand, cannot multistream, and might be more rotationally 
supported. Gas is also subject to shocks and galactic feedback, adding 
stochasticity, but also can tend to be smoother than the dark matter 
(Harford & Hamilton 201 1). 


An exotic implication of our result is in a superfluid Bose-Einstein 
condensate dark-matter scenario, with sufficient self-interactions that 
spin-driven dark-matter vortices (associated with spin on the scale 
of the whole halo) might form. Wave, or fuzzy, dark matter (VDM) 
(e.g., Hui 2021) is similar to this scenario, except that without strong 
self-interactions, all that can realistically form in y DM are small 
vortex loops of scale the de Broglie wavelength (Hui et al. 2021; 
Schobesberger et al. 2021). But an additional, topological obstacle 
for spin-driven vortices to form in a halo is that vortices cannot end 
in a point, but arise as loops. Any vortex associated with a halo's 
spin would have to complete the loop outside the halo. This becomes 
more plausible (but still seems far-fetched) if it can thread the cosmic 
web through connected filaments, associated with their spin. 


Filamentary gas since reionization is typically ionized, so if any 
processes might act differently on electrons and protons in filaments, 
e.g. in a Biermann battery mechanism (Kulsrud et al. 1997; Gnedin 
et al. 2000; Naoz & Narayan 2013), we speculate that rotation would 
generate a coherent magnetic field along the filament axis. A model 
with rotation in addition to shocks could help to understand the 
origin of cosmological seed magnetic fields. Also, an observational 
indication of a rotating filament could be a coherent magnetic field 
aligned with it, probed through e.g. synchrotron emission or Faraday 
rotation of a background polarized source (Brunetti & Vazza 2020). 


Our finding is consistent with the standard tidal-torque theory 
(TTT) of the origin and evolution of angular momentum in large- 
scale structure. Net angular momentum arises from velocities on the 
outskirts of collapsing structures, and gets transported to the centre 
by gravity. The net rotational velocity around a filament is compara- 
ble to the known quadrupolar pattern, revising the existing picture. 
Originating from the same large-scale environment, filaments and 
their nearby galaxies/haloes are expected to share a similar large- 
scale spin field, and this provides another reason for why the spin of 
galaxies might be correlated with that of their nearby filaments. 


We are unaware of a reason to think that our qualitative results, aris- 
ing from gravity and the TTT, might differ in any near-concordance 
ACDM cosmology; indeed, our Millennium and Illustris results are 
similar. But quantitatively, filament spins likely do depend somewhat 
on cosmological parameters. In the TTT, angular momentum grows 
with the linear-theory velocity field and the Universe's expansion, 
~ D3/2 (with D the growth factor) until ‘collapse,’ after which its 
angular momentum is conserved (White 1984). This holds also in 
2D, i.e. for infinitely long filaments (Neyrinck et al. 2020). Consider 
an ensemble of filaments, for which all properties are held constant 
except for the growth of velocities before filaments collapse. Their 
degree of spin would depend on the amplitude of the velocity field 
at collapse, i.e. on parameters like Oy and og. But the cosmic web 
and its filaments would also change with cosmological parameters, 
so the picture is likely not as simple as that, and worth further study. 


Our results indicate that the longest substantially rotating objects 
in the Universe are likely filaments. As we elaborate in $6, candidate 
*objects' for this title include cosmic-web components: dynamically- 
defined haloes, filaments, walls, and voids. Haloes, including the 
largest galaxy clusters, rotate, but filaments can be much longer. 
Walls and voids have similar lengths as filaments, and even substantial 
widths or depths. However, walls, and especially voids, are unlikely 
to rotate as substantially and coherently as the rotating filaments. 
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